Conduction through a quantum dot near a singlet-triplet transition 
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Kondo effect in the vicinity of a singlet-triplet transition in a vertical quantum dot is considered. 
This system is shown to map onto a special version of the two-impurity Kondo model. At any 
value of the control parameter, the system has a Fermi-liquid ground state. Explicit expressions for 
the linear conductance as a function of the control parameter and temperature T are obtained. At 
T — 0, the conductance reaches the unitary limit ~ 4e 2 //i at the triplet side of the transition, and 
decreases with the increasing distance to the transition at the singlet side. At finite temperature, 
the conductance exhibits a peak near the transition point. 
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The number of electrons TV in a Coulomb-blockaded 
dot is a well-defined quantity at low temperatures. The 
existence of a finite energy gap for excited states carrying 
a different charge, e(N ± 1), normally results in a sup- 
pressed low-temperature conductance through the dot. 
However, the suppression may be largely lifted Q| if N 
is odd. This phenomenon is explained in terms of the 
Kondo effect . In the case of tunneling through a spin- 
degenerate state, the conductance may reach the unit 
quantum 2e 2 /h, which corresponds to a perfect trans- 
mission H through the dot. The origin of the perfect 
transmission is in the formation of a "Kondo cloud" con- 
sisting of the itinerant electrons of the leads, which tends 
to screen the spin of the dot. The formation of such a 
collective state is accompanied by the appearance of the 
resonance for electrons right at the Fermi energy. Of 
course, this effect takes place only if the dot has a non- 
zero spin in the ground state. This is always the case for 
odd N. If N is even, then the spin of the dot is typically 
zero. However, if the spacing 5 between the last doubly- 
occupied and the following empty one-electron states in 
the dot is anomalously small, then the exchange inter- 
action, according to the Hund's rule, favors the triplet 
ground state of the dot Q . For a quantum dot formed 
in a GaAs heterostructure, 5 can be tuned by means of a 
magnetic field. Indeed, because of a very small electron 
effective mass, magnetic field applied perpendicular to 
the plane of electron gas has a strong orbital effect. At a 
certain critical value of 5 the singlet-triplet transition oc- 
curs. Tuning through such a transition was demonstrated 
recently in the experiments on vertical quantum dots [^) . 
Similar effect also occurs in lateral devices ||. The ap- 
plied magnetic field causes also spin splitting. However 
the Zeeman energy, which lifts the spin degeneracy, may 
remain sufficiently small, even compared with the Kondo 
temperature, because of a small value of the effective g- 
factor. Therefore the Kondo effect persists throughout 
the entire domain of parameters corresponding to the 
triplet spin state. 



The singlet-triplet transition was addressed in the re- 
cent work 0| , where quite special model of the quantum 
dot was considered. As it can be shown, the additional 
symmetries assumed in may result in a non-Fermi liq- 
uid behavior at the transition point. 

As we show below, a generic model of a quantum dot 
undergoing the singlet-triplet transition allows for a map- 
ping onto a the 2-impurity Kondo model. Using this 
mapping, we found that at finite temperature, the con- 
ductance G is a nonmonotonic function of the control 
parameter K with an asymmetric peak near the tran- 
sition point K = K* . The asymmetry becomes more 
pronounced at lower temperatures, and G(K) becomes 
constant (~ 4e 2 /h) at T = at the triplet side of the 
transition (K < K*), but decreases monotonously with 
K at K > K* . Despite the apparent non-analyticity of 
G(K), the system has a Fermi-liquid ground state for all 
K , including the transition point. 

The magnetic properties of a quantum dot can be un- 
derstood from the following Hamiltonian: 

ffdot = tndi s d ns - E S S 2 + E C (N- AO 2 . (1) 

ns 

Here, N = Y) ns d^ s (i„ s is the number of electrons in the 
dot, S = J2nss' d ns (~<?ss' /2) d ns i is the dot's spin, and the 
parameters Ec and E s are the charging and exchange en- 
ergies ||. We restrict our attention to the very middle 
of a Coulomb blockade valley with an even number of 
electrons in the dot (the dimensionless gate voltage AT 
is tuned to an even integer value). We assume that the 
spacing 5 between the last filled and first empty orbital 
states is of the order of E s , and that S is tunable, e.g., 
by means of a magnetic field B. In order to model the 
singlet-triplet transition in the ground state of the dot, 
it is sufficient to consider these two states; we will assign 
indices n = ±1 to them. The four low-energy states of 
the dot can be classified according to their spin S = 0, 1 
and it's z-projection S z . Labeling the states by these two 
quantum numbers, \S,S Z ), we find: 
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|ii) = < 1T rfl 1T |o) , |i-i) = 4 u d t _ u |o), 
|io) = -^ (<i/_ u + 4i At) 1°) . ( 2 ) 
|00) = dl lt dl u |0) 

where |0) is the ground state of the dot with N —2 elec- 
trons. Projected onto these states, the Hamiltonian of 
the dot becomes (up to a constant) 

VH dot V = ^ Yl l S ' SZ ) fei - 3 ^o) (5, 5 2 | , (3) 

s,s* 

where Kq = 8 — 2E S is the energy difference between the 
singlet and the triplet, and V is the projection operator 
on the system of states (|j). 

Upon the variation of magnetic field B, the singlet- 
triplet transition occurs at Kq = 0. Unlike the special 
case considered in Ref. @ , we are interested in the generic 
model with E s ^ at the transition point. 

In a vertical dot device, the potential creating lateral 
confinement of electrons most probably does not vary 
much over the thickness of the dot || . Therefore we as- 
sume that the electron orbital motion perpendicular to 
the axis of the device can be characterized by the same 
quantum number n inside the dot and in the leads; this 
quantum number is conserved in the process of tunneling. 
Thus, our model consists of the Hamiltonian of the dot, 
already discussed above, the Hamiltonian of the leads 
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and the tunneling Hamiltonian: 



(4) 



(5) 
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Here a = R,L for the right /left lead, and n = ±1 for 
the two orbitals participating in the singlet-triplet tran- 
sition; k labels states of the continuum spectrum in the 
leads, and s is the spin index. After a rotation in the R-L 
space, 
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the 4> field decouples: H T = J2 n ks tni>l, ks dns + H.c; here 
t 2 n = t\ n + t 2 Rn . Next we integrate out the virtual transi- 
tions to the states with M ± 1 electrons by means of the 
Schrieffer- Wolff transformation. The resulting effective 
low-energy Hamiltonian includes the operators 



-d n i a >V. 



The effective Hamiltonian may be conveniently written 
in terms of two fictitious 1/2-spins Si, 2 which represent 



the same symmetries as the set of states (Q). This one- 
to-one correspondence between the basis states allows us 
to represent operators S nn / in terms of S 12 - We find the 
following relations: 



_( Sl -s 2 ) = 7 =s_, 



£ 



^ ^ S— n,? 

n 

S_n,n = V2[Si X S 2 ] = V2T 



(6) 



Using (^J) , the effective Hamiltonian is written in a form, 
resembling M the two-impurity Kondo model Gof : 



H = ^2^l ks Tpnks+K(S 1 -S 2 ) 

nks 



H n = J„ (s„„ ■ S+) + Vnp nn (Si • S 2 ) 



(7) 
(8) 



+-j= [(s_„,„ • S_) + 2m (s_ n ,„ • T)] . 



Here the particle and spin densities in the continuum are 

Ifs 



Pnn = £ ^Ls^k' 



ns-f ^nn' 



£ 

kk' ss' 
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and the bare values of the coupling constants are 

7 - 1 - 2t + it - 1 v - t2 +i + t -i (q) 

We did not include in Eq. (m) the Hamiltonian of the 
</> field, and other terms which are irrelevant for the low- 
energy renormalization. The Schrieffer- Wolff transfor- 
mation also produces corrections to Ko, so K does not 
coincide with its bare value. This difference is not impor- 
tant, as it only affects the position of the singlet-triplet 
transition, but not its nature. A common factor / in the 
last two terms of Eq. (|§|) comes from the conservation of 
the orbital index n, see Eq. (||). 

To simplify the analysis of Eqs. (@)-@, we further re- 
strict our attention to the symmetric case t an = t a , for 
which the definition (|J) reduces to 



J n = J = I = 2V = 2 (<| + t 2 R ) /E c 



(10) 



This simplification is adequate to the experimentally rel- 
evant case of very thin barriers separating the dot from 
the leads, and, more importantly, only insignificantly af- 
fects the low-energy properties of the model. To calcu- 
late the differential conductance in the leading logarith- 
mic approximation, we apply the "poor man's" scaling 
technique . The procedure consists of a successive in- 
tegration out of the high-energy degrees of freedom, and 
yields the set of scaling equations 



dJ 1 T 2 



on dl , . dV o 

I*),- = 2vI(J + V), 1E =2vI>, 



(11) 
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with the initial conditions (|To|). Here £ = hiEc/D, and 
v is the density of states in the leads; the initial value of 
the energy cutoff D is D = Ec ■ This procedure also gen- 
erates non-logarithmic corrections to K. In the following 
we absorb these corrections in the re-defined value of K. 
Equations (0) are valid for D > \K\,T. 

At certain value of £ = Co, the inverse coupling con- 
stants simultaneously reach zero: 



l/J(£ ) = l/I(£ ) = l/V(£ ) = 0. 

This defines, through the equation £q = hiEc/To, 
characteristic energy scale of the problem: 



the 



T = E c cxp [-t /vJ{0)}. 

Here To is a parameter that depends on the initial con- 
ditions and should be found numerically. We obtained 
t = 0.36 (see Fig.|l|). Thus, the strong coupling regime 
is reached at much higher temperatures, than in the usual 
Kondo model (for which the Kondo temperature would 
be given by the same expression as To, but with to = !)■ 




FIG. 1. Numerical solution of the scaling equations. The 
RG equations ( [ill ) are rewritten in terms of the new variable 
r = vj(0) InEc/D and the new functions x(t) = J(0)/J(t), 
y(r) = 7(0)/7(r), z(r) = V(0)/V(t). The three functions 
reach zero simultaneously at r = to = 0.36. 

The solution of the renormalization group (RG) equa- 
tions dll] ) can now be expanded near £ — £q- To the 
first order in £o — £ = InD/To, we obtain 

vJ vi 2vV 

where A = 2 + a/5 ~ 4.2 is a model-independent constant, 
i.e., it does not change if the restriction t an = t a is lifted. 

The solution ( jl2| ) can be used to calculate the differ- 
ential conductance at high temperature T ^> \K | , Tq. In 
this regime, the coupling constants are still small, and 
the conductance is obtained by applying a perturbation 
theory to the Hamiltonian ((?])- (||) with renormalized pa- 
rameters. This yields 



G/G =A\n- 2 T/T , 



(13) 



A= (3^/8) (A + ir 2 1 + A+(A-1) 
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is a numerical constant, and 



Go 



4e2 
~h 



2t T .t 



ft + tl 



(14) 



Note that ( |13|) includes contributions from both the pro- 
cesses conserving the orbital index [the first two terms 
in Eq. (0)], and the processes involving an inter-orbital 
scattering. 

Away from the singlet-triplet degeneracy point, \K\ > 
T , the RG flow yielding Eq. (|l|) terminates at energy 
D ~ \K\. At the triplet side of the transition (K < 0), 
the two spins Si,2 are locked into a triplet state. The sys- 
tem is described by the effective 2-channel Kondo model 
with 5 = 1 impurity, obtained from Eqs. (0)-(j|) by pro- 
jecting out the singlet state: 



-^triplet = £,ktpi ks 1pnks + Jt*^2 ( Sn ™ ' S) + Vi 



np rl 



nks 



here J t = J (\K\) , V t = V (\K\) /A. As D is further low- 
ered, Jt is governed by the standard Kondo RG equation 

dJ t /d£ = vJ?, £ = InD/ \K\ . 

The solution of this equation, l/v.J t (D) — In D/T^, is 
expressed through the independent Kondo temperature 
Tc (K) — \K\ exp (— l/vJ t ), which, using (|T2]), is in turn 
expressed through Tq as 



T fe /To = (T / \K\Y 



(15) 



Recall that the exponent here, A « 4.2, is universal. 
Eq. (|l^) was obtained also in 0, but with a rather dif- 
ferent value of the exponent A (according to 0, A < 1 
and appears to be non- universal) . 

Since Vt is not renormalized, the differential conduc- 
tance at T <C — K is dominated by the exchange (the 
second term in iTtripiet) and is given by 



G/G = f (T/T fc ) = / 



T 
To 



\K\ 
To 



(16) 



where / (x) is a smooth function that interpolates 
between / (x > 1) = (7r 2 /2) In" 2 x and / (0) = 1. 
It coincides with the scaled resistivity /(T/Tr-) = 
p(T/Tfc)/p(0) for the symmetric two channel S = 1 
Kondo model. The conductance at T = (the unitary 
limit value), Go, is given above in Eq. (|i~4|), see also Fig. g. 

On the singlet side of the transition, K > T , the scal- 
ing terminates at D ~ K, and the low-energy effective 
Hamiltonian is 



where 



H. 



singlet 



V s ^2 np r , 



3 



where V s — — 3V (K) /4. The conductance at tempera- 
tures T <C K is given by 



a/G, = uhr 2 K/T , r, ( i^^-i 



0.5. (17) 



only at one special value of magnetic field, such that the 
Zeeman energy equals S. Accordingly, G (T, B) exhibits 
a Kondo peak at this strength of the field at all temper- 
atures. 



The only regime of parameters left beyond the range 
of the validity of the equations (13), and (0) is 

that of the very vicinity of the singlet-triplet transition 
K = at low temperature T < Tq. Understanding of 
this regime requires knowledge of the properties of the 
system's ground state. This can be inferred from the fol- 
lowing simple symmetry-based argument, devised origi- 
nally for the two impurity Kondo model |l(| . The Hamil- 
tonian (^), (|^) in the symmetric case J n — J is invariant 
with respect to the particle-hole transformation 

ipn,k,s -> sip 1 _ n _ k _ s . (18) 

Consider now the scattering problem at the Fermi energy: 

C=W^^- (19) 

Particle-hole symmetry ( |l8| ) implies that for any $ ln ( out ) 
which solves Eq.@, fc^W) = i£$out(in) wiU be a g0 _ 
lution as well (here E nSi „/ s / = ^L' "™') an d cr l are the 
Pauli matrices). Substituting this into ( |l9| ) one finds a 
relation for the S-matrix: (XS) 2 = 1. Among the diag- 
onal S, this is satisfied by S nSi „' S ' = 6 nn /5 ss >e~ 2me with 
arbitrary 9. In other words, the particle-hole symmetry 
imposes a restriction on the scattering phase shifts at 
the Fermi energy: 9 ns — nO. The phase 9, and therefore 
scattering phases 9 ns vary continuously from ±7r/2 to 
as K is varied from K — — oo to K = +oo. At the triplet 
side of the transition, — K To, the low-temperature 
physics is described by the effective 2-channel S = 1 
Kondo model, as discussed above. In that model, the 
zero-energy phase shift 8 = tt/2. This means that there 
exists K* > -T such that 9(K) = tt/2 for all K < K* . 
On the other hand, at the singlet side of the transition, 
0(K) decreases with K logarithmically at large positive 
K ^> To, but approaches tt/2 at smaller K because of 
larger values of V s in -ff s ingict- We expect therefore that 
9{K) is a continuous, though a non-analytical, function. 
The singlet-triplet transition may be associated with the 
point of the non-analyticity K = K* . According to these 
arguments, it is natural to conjecture that the Fermi liq- 
uid description jlj] is applicable at all K and thus the 
conductance at T = is G = Go sin 2 9{K). The inter- 
orbital scattering processes make no contribution to G in 
this regime, unlike at high temperature T 3> | if | , To . 

It should be noticed that the problem under considera- 
tion is quite different from what one faces if the Zeeman 
splitting is the leading effect of the magnetic field |T5| ] 
(which may happen if magnetic field is applied in the 
plane of the lateral dot). In this case, the two- fold de- 
generacy of the ground state of an isolated dot appears 



-To ° ^0 T K 

FIG. 2. Tf-dependence of the differential conductance at 
two different temperatures. The two asymptotes merge at 
K > T,T . 

To conclude, we considered the linear conductance 
G(T, K) of a vertical quantum dot near the singlet-triplet 
transition. The transition occurs due to the strong or- 
bital effect of an external magnetic field. At high tem- 
perature G exhibits a peak near the transition point, in 
agreement with the experiments |5| . At low temperature 
G reaches the unitary limit value at the triplet side of 
the transition, and decreases monotonously at the sin- 
glet side (see Fig. ^). The characteristic energy scale of 
the effect is much larger than that for the usual Kondo 
effect with the similar system parameters. 
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